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Abstract 

The refractive index of a dielectric medium comprising both passive and inverted components 
in its permittivity was determined using two methods: (i) in the time domain, a finite-difference 
algorithm to compute the frequency-domain reflectance from reflection data for a pulsed plane 
wave that is normally incident on a dielectric half-space, and (ii) in the frequency domain, the 
deflection of an obliquely incident Gaussian beam on transmission through a dielectric slab. The 
dielectric medium was found to be an active medium with a negative real part for its refractive 
index. Thereby, a recent controversy in the scientific literature was resolved. 

Keywords: Active medium; Negative refraction; Time-domain analysis; Frequency-domain anal- 
ysis 

1 Introduction 

A complex number z = ip e^'^ possesses two square roots: ^/z = ztly^l e*^"^/^^ . This presents a 
problem when a physical quantity (e.g., refractive index) is expressed as the square root of another 
physical quantity (relative permittivity). Which root is the physical one? 

Our work grows out of a recent controversy over this question in the context of negative re- 
fraction of electromagnetic plane waves. Chen, Fischer, and Wise (CFW) considered an isotropic 
nonmagnetic medium with a relative permittivity scalar 
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e (w) = 1 ^ aeujj [uJj - {u; + iPiujif] , (1) 
e=i 

as a function of the angular frequency lv [1]. Herein, the constants ai = 2.4401, a2 = —0.14348, 
Pi = 0.028571, P2 = 0.020000, oji = 2.6371 x lO^^ rads'^ and uj2 = 3.7673 x lO^^ rads'^. 
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The permittivity scalar comprises a passive and an inverted component. Most importantly, for 
Ao G [445, 535] nm, where Aq denotes the free-space wavelength, Im(e) < 0, as shown in Fig. 1; here 
and hereafter, an exp(— iwt) time-dependence is implicit for all frequency-domain field phasors. 
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Figure 1: Real (solid) and imaginary (dashed) parts of e plotted against free-space wavelength (in 
nm). 

CFW deduced that for Aq G [445, 535] nm, the refractive index n = \f\ must be such that 
Re(n) < and Im(n) > 0. This deduction stemmed from the premise that the phase angle of e is a 
continuous function of Aq, which has been supported by certain recent theoretical arguments [2, 3]. 
Accordingly, CFW argued that their medium could refract negatively. 

Alternatively, on the basis that (i) energy flow grows along the direction of propagation when 
Im(e) < [4]; or (ii) the real part of the wave impedance is positive- valued on the grounds 
of causality [5], it may be deduced that Re(n) > and Im(n) < for the CFW medium for 
Ao G [445,535] nm [6, 7]. CFW disputed this alternative view [8, 9]. 

In order to resolve the issue, two studies were undertaken. First, following a procedure adopted 
by Wang and Lakhtakia [11], we performed a time-domain calculation that did not explicitly invoke 
the frequency-domain concept of refractive index. We used a finite-difference algorithm to solve 
the time-domain Maxwell equations for a pulsed plane wave reflected, at normal incidence, from a 
half-space filled with the CFW medium. Then we transformed the time-domain electric field of the 
reflected pulse to the frequency domain and computed the reflectance as a function of Aq. Second, 
in the frequency domain, we considered a Gaussian beam propagating through a slab of the CFW 
medium at an oblique angle. The reflection and transmission coefficients were computed — without 
utilizing the refractive index — by solving the reflection-transmission problem as a boundary- 
value problem, and the deflection of the transmitted beam with respect to the incident beam was 
determined. 
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2 Time— domain analysis 



Let us begin with the time-domain study. Suppose the CFW medium occupies the half-space 
z > zl, {zl > 0), and possesses the time-domain relative permittivity 



e{t) = S (t) + ^ a^uji exp {-fiiUft) sin {uet) U (t) , 
e=i 

where 6{t) is the Dirac delta function, U{t) is the unit step function; note that 

/•CO 

e [oo) = e{t) exp (iojt) dt . 



(2) 



(3) 



The other half-space z < zl is vacuous. 

A pulsed plane wave propagating in the +z direction is introduced at z = 0, so that for t > 



E = 0, = X 



exp 



\V2toJ 



cos (coct) . 



(4) 



The electric field E{z,t) = xE{z,t) is polarized along the x axis, the magnetic field H(z,i) = 
yH{z,t) is polarized along the y axis, 770 = ^/no/eo is the intrinsic impedance of free space 
(permittivity eo and permeability /iq), cq = 1/^eo/xo is the speed of light in free space, Ut sets the 
energy density of the pulse, tq is the time constant, td the time delay, and Uc the carrier frequency. 

Upon writing the components of the electromagnetic field in a column 2-vector as [F_{z,t)] = 
[E {z, t) , H {z, t)]^, where the superscript indicates the transpose, and substituting the foregoing 
expressions for e {t) into the Maxwell curl postulates, we found the matrix partial differential 
equation 



dt [F{z,t)]=co[Z]d,[F {z, t)] - eo ' dt [ Q {z, t)] 
for z,t > 0. In this equation, 

\V] = \ ° '"f" 

is the vacuum propagation matrix, the column vector 

[Q{z,t)]=eo f[W{z,t')] [F{z,t-t')]dt', 
Jo 

and the matrix [lK(2:, i)] is null-valued for z < 2;^,, but 



(5) 
(6) 

(7) 



(8) 



for z > Zl- The upper limit on the right side of eqn. (7) accounts for [F_{z,t)\ being null-valued 
for t < 0. 

We computed the spatiotemporal evolution of the pulsed plane wave over the domain {{z,t)\z ^ 
[0, zr] , z/j > ZL-,t > 0}, which was discretized into space steps of length A2; and time steps of 
duration At = /3Az/co, where /? < 1 is a stability parameter. We discretized eqn. (5) and solved 
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Figure 2: Real (dashed) and imaginary (dotted) parts of the indexes of refraction, as defined in the 
frequency domain. 



it using over the chosen domain using a leapfrog finite-difference algorithm. Further details of our 
solution procedure are reported elsewhere [10]. 

We chose zq = fj,m, zl = 20 /xm, and zr = 40 ^m. The electric field of the incident pulse 
Ei (t) was recorded at z = /um, and that of the reflected pulse Er (t) at 2; = 18 fim. The incident 
pulse had the parameters C/f = 1 J m^^, tq = 1 fs, = 5 fs, and lOc = 3.8838 x 10^^ rads~^, 
so that its bandwidth was centered at the free-space wavelength Aq = 485 nm with a full-width 
half-maximum of about 220 nm. Then, we used the fast Fourier transform (FFT) to approximate 
the corresponding phasors E^ (Aq) and Er (Aq) and found the reflectance from the time-domain 
calculations as 



Rt (Aq 



Er (Aq) Ei (Aq) 



(9) 



There are two possibilities for the refractive index n, viz., n± such that n_ = — and Im(n+) > 
0. Plots of Re {n±) and Im {n±) as functions of Aq are shown in Fig. 2. From these two possibilities, 
we found the frequency-domain reflectances 
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i?±(Ao)= [n±(Ao)-l][n±(Ao) + l]- 



(10) 



for plane waves normally incident on a half-space. As = 1, Rt can be used to distinguish 

between them if R± 7^ 1 [11]. 

Plots of Rt, and R- vs. Aq are shown in Fig. 3 over the bandwidth covered by the incident 
pulse. As l-Rtj > 1 over at least part of that bandwidth, the CFW medium is active over that 
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part of that bandwidth. However, we found that neither the pulse reflected from nor the pulse 
refracted into the medium grew unboundedly. Furthermore, as the reflectance Rt obtained from 
the time-domain calculations closely matches the refractive index with positive imaginary part 
(i.e., n+) is the correct one. 
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Figure 3: Reflectances Rt (blue solid), R+ (red dotted), and (green dashed). 



3 Frequency— domain analysis 



Next we turn to the frequency-domain study which does not require the specification of n as either 
n_ or n-i-. Suppose that the half-space is now replaced by a slab of thickness L, as schematically 
illustrated in Fig. 4. The slab — which consists of dielectric material with relative permittivity 
defined in eqn. (1) — is sandwiched by two vacuous half-spaces. 
A 2D beam with electric field phasor [12] 



Ej {x,z,Xo) 



3i(??)^(??) e^'^+('')-'" di}, 



(11) 



for z < 0, is incident upon the slab at a mean angle 9i relative to the slab normal direction z. The 
beam is represented as an angular spectrum of plane waves, with 



k±(i?) 



kn 



^ COS Oi + \/l -i92 



sm I 



cos ( 



(12) 



where fco = 2tt/\q. The angular-spectral function ^{-d) is taken to have the Gaussian form [12] 

(13) 



y^w) = — j= exp 



27r 



1 2 
-- (/co Wq l}) 



with wq being the width of the beam waist. Two polarization states are considered: parallel to the 
plane of incidence, i.e.. 



ei(i9) 



x 



■d sin 9i — \/l — cos 6, 



+ (^'d cos 9i + \/l -i?2 

sin 9j I 



(14) 
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Figure 4: A beam is incident onto a slab at a mean angle Qi with respect to the unit vector z 
normal to the planar interface. The incident beam strikes the slab at the coordinate origin (i.e., 
a; = 0, y = and z = 0). The transmitted beam emerges from the slab at z = L at a point with (a) 
X > if the refractive index of the slab is positive (dashed red arrow); or (b) x < if the refractive 
index of the slab is negative (solid red arrow) . 



and perpendicular to the plane of incidence, i.e., 

ei(??) = = y. 



(15) 



As the incident beam has the spatial Fourier representation (11), the reflected and the transmit- 
ted beams must also have similar representations. The electric field phasor of the reflected beam 
is given as 



(x,2;, Ao) 



for z < 0, with 



er(??) = < 



+ (i?cos6'i + Vl-i?2 sin Qi 
for ej('!9) = e|| 
e_|_ for ej(t?) = e_|_ 
The electric field phasor of the transmitted beam is given as 



?? sin Bi - Vl -192 



cos f ,• X 



Ej (x,2;, Ac 



(i?)^(i?)e*^+('')-('-^^) 



et(i?) 



(16) 



(17) 



(18) 



(19) 



for z > with 

' t|| e0) for ei{'&) = e,, 
t_L e_L for ei{-d) = 

The reflection coefficients ry _|_ and transmission coefficients ty^^ were calculated as functions of "d 
by solving a boundary-value problem [13], as described in the Appendix. 
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We fixed the mean angle of incidence 9i = 60°, the free-space wavelength Aq = 485 nm, the 
beam waist wq = 1.75Ao, and the slab thickness L = 4Ao. The restriction G [— 1, 1] was imposed 
to exclude evanescence. The numerical values for the beam waist and slab thickness were chosen 
in order to accentuate the clarity of Fig. 5, which shows the energy density in both half-spaces, as 
defined by 



for z/Xo G (-8, 12) and x/Xq G (-25, 25). 

As illustrated in Fig. 1, at Aq = 485 nm the relative permittivity of the CFW material is 
e = 0.51 — 0.87i. The corresponding reflection and transmission coefficients were numerically 
determined at i? = as ry = -3.16 + 0.58z, r± = -1.14 - l.Sli, t\\ = (3.69 + 3.66z) x IQ-^ 
and t± = (—1.82 — 0.10) x 10~^. In order to make visible the tiny fraction of the beam that is 
transmitted, the values of |Etp in Fig. 5 have been amplified by a factor of 5 x lO"*^^ for ej = ey, 
and by a factor of 3 x 10^^ for ei = e±. The fact that the CFW medium is active at Aq = 485 nm 
is clear from |r||p > 1 and |r_Lp > 1 (in fact, |r||p = 10.32 and |r_Lp = 3.00). 

From Fig. 5 we conclude that the beam undergoes negative refraction at the two interfaces 
between the CFW slab and free space. While this occurs for both polarization states, it is more 
noticeable for the parallel polarization state. 

4 Concluding remarks 

In conclusion, there is an ambiguity inherent in the frequency-domain concept of refractive index, 
concerning the choice of square root. This pertains to both the bending of light at planar interfaces 
and the determination of whether the medium under consideration is active or passive. Using (i) 
a time-domain method and (ii) a frequency-domain method, neither of which explicitly invokes 
the refractive index, we resolved these issues for the two-component CFW medium characterized 
by the relative permittivity given in eqn. (1). The CFW medium was found to be (a) an active 
medium — contrary to the claims of CFW [1]; and (b) negatively refracting — contrary to other 
recent claims [6, 7]. 

Acknowledgements: We thank S. A. Ramakrishna for helpful discussions. JBG is supported by 
a Beckman Postdoctoral Fellowship. TGM is supported by a Royal Society of Edinburgh/Scottish 
Executive Support Research Fellowship. 



Appendix 

The reflection coefficients r[[^^ and transmission coefficients t\\^± are straightforwardly calculated 
by solving a boundary-value problem as follows [13]. 




for z > L 



for z < 



(20) 
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Figure 5: Normalized |Ep is mapped in the xz plane for a 2D Gaussian beam incident onto a CFW 
dielectric slab at a mean angle 9i = 60°. Ej is polarized parallel (top) and perpendicular (bottom) 
to the plane of incidence. The red line indicates the mean beam position in the absence of the 
dielectric slab. 



Consider the plane wave with electric and magnetic field phasors 

E(x,z) =e{z,6) eyLp (ik^xsinO) 
H(x,z) ='h{z,9) exp (ifcoxsin^) 

propagating in the xz plane and incident on a dielectric slab with relative permittivity e occupying 
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the region between z = and z = L. The angle 9 is related to the mean angle 9i of a Gaussian 
beam and the parameter by the twin relations 

sin = •& cos Oi + Vl — "f?^ sin 9i | 



cos 



e = -'& smOi + Vl^^ cos d^,- 



(22) 



Substitution of eqn. (21) into the source-free Maxwell curl postulates VxE(a;, z)—i(jJiiQil{x, z) = 
and V x H(x, 2:)+iu;eoeE(x, z) = delivers four differential equations and two algebraic equations. 
The latter two equations are easily solved for and hz- Thereby, the four differential equations 
may be expressed in matrix form as 

d 



dz 



where 



is a column vector and 



l{z,e) = ex{z,e), ey{z,e),hx{z,e),hy{z,e) 



m 



is a 4x4 matrix with 












mp 








-m 





























sin^ 9 






p = l 


e 







The solution to eqn. (23) is conveniently expressed as 



f{L,9) 

in terms of the transfer matrix [13] 



M(L,0) /(0,e) 



MiL,9)\=Y, 



£\ 



iko 



m\ ^} • 



(23) 
(24) 

(25) 

(26) 
(27) 

(28) 



,z < 0. 



(29) 



Since the evaluation of M_{L-, d) ^ ^ power series does not invoke the refractive index, ambiguities 
associated with the determining the correct square root of e are avoided. 

Now we turn to the incident, reflected and transmitted plane waves. Let the incident plane 
wave be represented in terms of linear polarization components as 

ei{z, 9) = [a_L y + a\\ (sin 6* z — cos 9 x)] 
X exp {ik^z cos 9) 

h.i{z, 9) = [a_L (sin^z — cosmic) — oy y ] 
xr/^^ exp {ik^z cos 9) 

The corresponding reflected and transmitted plane waves are given as 

er{z, 9) = \a±r± y + ayrn (cos^x + sin^z)] 

X exp {—ikoz cos 9) 

hr{z,9) = [a±r± (cos x + sin ^ z) — a||r|| y] 
xrjQ^ exp (—ikoz cos 9) 



,z<0 



(30) 
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and 



et{z, 0) = \_a_\_t_\_ y + a||t|| (306* z — cos 9 x)] 
X exp [iko{z — L) cos 0\ 

h.t{zj 0) = \a±t± (sin^z — cos^x) — a\]t<^\ y ] 
y^rjQ^ exp [ikQ{z — L) cos 6] 



,z>L, 



(31) 



respectively. By application of the boundary conditions at z = and z = L to the solution (27), 
the reflection and transmission coefficients are found to be related by the matrix algebraic equation 



[Ki9)][t±,t\\,0,0f = M{L,e) [Ki9)][l,l,r±, 



(32) 



wherein 








— cos 9 





cos 9 




im)] = 


1 

—rjQ^ cos 6 






1 

tiq^ cos 9 







(33) 












-^0-^ . 





Thus, the reflection and transmission coefficients emerge as components of the 4x4 matrix 



S 



[Kie)]'' Km 9) [Km^ 



as per 



r_L = - 


S 


3.( 


S 










s 


42 ( 


s 











44 



-1 



-1 



t± = 

til = 
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11 
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13 
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22 


s 


24 


s 
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s 



33 > 



(34) 

(35) 
(36) 
(37) 
(38) 
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